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The current study proposes a model of the cardiovascular system that couples heart cell 
mechanics with arterial hemodynamics to examine the physiological role of arterial blood 
pressure (BP) in left ventricular hypertrophy (LVH). We developed a comprehensive 
multiphysicsand multiscale cardiovascular model of the cardiovascular system that 
simulates physiological events, from membrane excitation and the contraction of a cardiac 
cell to heart mechanics and arterial blood hemodynamics. Using this model, we delineated 
the relationship between arterial BP or pulse wave velocity and LVH. Computed results 
were compared with existing clinical and experimental observations. To investigate the 
relationship between arterial hemodynamics and LVH, we performed a parametric study 
based on arterial wall stiffness, which was obtained in the model. Peak cellular stress of the 
left ventricle and systolic blood pressure (SBP) in the brachial and central arteries also 
increased; however, further increases were limited for higher arterial stiffness values. 
Interestingly, when we doubled the value of arterial stiffness from the baseline value, the 
percentage increase of SBP in the central artery was about 6.7% whereas that of the 
brachial artery was about 3.4%. It is suggested that SBP in the central artery is more 
critical for predicting LVH as compared with other blood pressure measurements. 

Keywords: Computer Simulation; Blood Pressure; Hypertrophy, Left Ventricular; 
Integrative Cell-System-Arterial Network Model 



INTRODUCTION 

There have been a number of recent studies examining the re- 
lationship between arterial wall stiffening and cardiac cellular 
overload (1, 2). However, despite substantial progress in under- 
standing this issue, the detailed mechanisms underlying this 
relationship remain unclear. Therefore, a novel approach is re- 
quired to delineate the causal relationship between arterial 
stiffness and heart hypertrophy. The use of computer simula- 
tions with mathematical models is an alternative method and 
has been widely used for the analysis of cardiovascular system 
dynamics (3-6). 

Lumped models of the cardiovascular system are very simple 
and can be easily coupled with the nervous system for long- 
term simulation (7). However, the model cannot accurately 
predict microscopic phenomena in cells or tissues. To over- 
come this limitation, we proposed a cell-system integrated 
model of the cardiovascular system that included cardiac cells. 



the Laplace heart, and the vascular system (8-10). However, the 
cell-system model uses relatively a simple lumped model for 
the arterial system and lacks detailed information regarding ar- 
terial pressures and pulse wave velocity (PWV). On the other 
hand, mathematical models focusing on arterial hemodynam- 
ics have been proposed by several groups (11-13). These mod- 
els have the ability to represent the spatial distribution of arteri- 
al hemodynamics and PWV However, the time-varying capaci- 
tance heart model (14) was incorporated into these models by 
which the relationship between heart cellular mechanics and 
arterial hemodynamics could not be analyzed. 

In this study, the arterial network model of Ozawa et al. (11) 
was incorporated into our previous cell-system model (8) to 
delineate the physiological relationship between arterial pres- 
sure and left ventricular hypertrophy (LVH). To verify the pres- 
ent method, we calculated the effects of arterial stiffness, arteri- 
al hemodynamics, and pulse wave velocity on cardiac cellular 
mechanics, and compared the results with experimental and 
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clinical observations. Using this methodology, the physiologi- 
cal effects of arterial stiffness variation on LVH were assessed in 
an integrated manner. 

MATERIALS AND METHODS 

To simulate heart mechanics and arterial pulse waves at a mul- 
tiscale level, we developed an integrative mathematical model 
that combined cell excitation- contraction coupling with heart 
mechanics, system circulation, and hemodynamics of the arte- 
rial network. In our previous studies (8, 10), we combined the 
ventricular cell model with the Negroni and Lascano model (NL 
model) (15) and eventually with system circulation. Here, we 
used the same approach for the atrial and ventricular models, 
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Fig. 1. Schematic of the integrative cell-system simulation model designed for effec- 
tive multiscale numerical analysis from cellular to tissue mechanics. BLH, biological 
Laplace heart. 



but with the addition of arterial network branches (Fig. 1). 

As described in our previous report (10), the present compu- 
tational approach to cardiac cells is based on the electrophysio- 
logical models of human atrial and ventricular myocytes pro- 
posed by Nygren et al. (16) and ten Tusscher et al. (17) (TNNP 
model), respectively. In addition, these cellular electrophysio- 
logical models were combined with the NL model of cross- 
bridge dynamics (CBD) for contraction (18). Cross-bridge dy- 
namics in the NL model are described by a four-state system 
consisting of free troponin (T), Ca^^-bound troponin (TCa), 
Ca^^-bound troponin with attached cross-bridges (TCa*), and 
troponin with attached cross-bridges (T*). The total force (F) of 
a single muscle unit is given by the following relationship: 

F = Fb-HFp (a) 
where Fb is the equivalent cross-bridge force generated by the 
contribution of all cross-bridges attached in parallel within the 
muscle unit and Fp is the force developed by the elastic element 
of the muscle unit. The cross-bridge force generated in (a) is 
presented as follows, according to the NL paper (15): 
Fb = A-([TCa1-H[Tl)-h (b) 
where A is the model parameter and h is the cross-bridge elon- 
gation length. The related parameters in the NL models are 
presented in Table 1. 

The forces obtained from the cell-NL models were converted 
into atrial and ventricular pressures by applying Laplace's law 
to the atrial and ventricular models. In our previous report, we 
referred to this type of heart model as a biological Laplace heart 
(BLH) (9). The NL, as well as the spherical atrial and ventricle 
models, were explained in detail in our previous reports (8, 9). 

We used a lumped parameter model to simulate a cardiovas- 
cular system. With the exception of the heart and arterial com- 
partments, the hemodynamic elements were the same as those 
described by Heldt et al. (19). The entire model is shown in Fig. 



Table 1. Equations and parameter values for the cross-bridge dynamics model 



Heart chamber 


Forces 


Equations & parameter values 


References 


Left ventricle 


CBCF 

Cellular passive elastic reaction force 


Fb = A([TCa1-F[T1)h 

A = 1,800 mN/mm7|jm/mM 

([TCa*]: concentration of cross-bridge attached to calcium-bound troponin, 
[T*]: concentration of cross-bridge attached to troponin, 
h: cross-bridge elongation length (CBEL)) 

Fp = E(e°(L^L^-^>-l), (L>Le) 

(L: half sarcomere length, Le: half sarcomere length at equilibrium state) 
E = 2pM/s, D = 10 
Fp = -B(1-L/Le),(L<Le) 
B = 20 


Negroni and Lascano (1 5) 
Negroni and Lascano (1 5) 

Landesberg and Sideman (29) 
Landesberg and Sideman (29) 


Left atrium 


CBCF 

Cellular passive elastic reaction force 


Fb = A([TCa1-Fri)h 

A = 450 mN/mm7|jm/mM 

Atrial cell -generated maximum force is 1/4 of that of ventricular cell 

Fp = E(e°(^L^-^>-1), (L>Le) 

E = 1 .0 |jM/s, D = 20 

Atrial cells are stiffer than ventricular cells 

Fp = -B(1-LyLe),(L<Le), 

B = 20 


Negroni and Lascano (15) 
Beyar and Sideman (30) 

Landesberg and Sideman (29) 
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Venous circulation model 
Upper body 



Arterial circulation model 




Fig. 2. Schematic of the integrative vascular system model, up, upper body; sup, superior vena cava; pa, pulmonary artery; pv, pulmonary vein; ro, right ventricular outflow; la, 
left atrium; Iv, left ventricle; inf, inferior vena cava; ab, abdominal vena cava; kid, kidney; sp, splanchnic; II, lower limbs; th, thoracic; rv, right ventricle. 



2. The peripheral circulation is divided into the upper body, re- 
nal, splanchnic, and lower extremity sections; the intrathoracic 
superior, inferior vena cava, and extrathoracic vena cava are 
identified separately. The model thus consists of compart- 
ments, each of which is represented by a linear resistance (R) 
and compliance (C) that can be either linear, nonlinear, or time- 
varying. The NL model and Laplace's law were applied to calcu- 
late heart mechanics. All of the related parameters in the sys- 
temic circulation model are presented in Table 2. 

The arterial network model proposed by Ozawa et al. (11) 
was incorporated into the current cell-system model system as 
shown in Fig. 2. Model of Ozawa et al. was based on the numer- 
ical solution of one-dimensional equations for fluid and vessel 



wall motion in a geometrically accurate branching network of 
the arterial system (28 elements), including energy losses at bi- 
furcations. Blood flow and pulse wave propagation in each ele- 
ment were described by one-dimensional Navier-Stokes equa- 
tions with axial velocity, fluid pressure, and vessel cross-sec- 
tional area as variables (14). These equations were solved using 
the finite difference method on one-dimensional meshes for 
each element of the arterial network branches (14). Rebound of 
arterial waves from terminal or bifurcation branches were nu- 
merically implemented through boundary conditions. As simi- 
lar to Ozawa et al. (11), we used the method of charateristics for 
the boundary conditions. Detailed specifications of the four 
large arteries in the arterial network are shown in Table 3. The 
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Table 2. Lumped parameter values for the cardiovascular system model 



Locations 


Elements 


Parameter values 


References 


Left atrium 


Outflow resistance 


Rav = 0.0025 PRU 


Ursino and Innocenti (31) 






( = mmHgs/mL) 






Inflow resistance 


Rpv = 0.006 PRU 


Ursino and Innocenti (31) 


Left ventricle 


Outflow resistance 


Rio = 0.004 PRU 




Aorta 


Outflow resistance 


Rio = 0.03 PRU 






Compliance 


0.03 mL/mmHg 




Peripheral 


Rupi 


3.9 PRU 




resistance 


Rup2 


0.23 PRU 






Rsup 


0.06 PRU 






Rpv 


0.01 PRU 






Rro 


0.003 PRU 


Heldt et al.(19) 




Rinf 


0.015 PRU 






Rab 


0.01 PRU 






Rkidl 


4.1 PRU 






Rkicl2 


0.3 PRU 






Rspl 


3.0 PRU 






Rsp2 


0.18 PRU 






Rin 


3.6 PRU 






R|I2 


0.3 PRU 





PRU, peripheral resistance unit (mmHg s/mL); up, upper body; sup, superior vena 
cava; pv, pulmonary vein; ro, right ventricular outflow; inf, inferior vena cava; ab, ab- 
dominal vena cava; kid, kidney; sp, splanchnic; II, lower limbs. 



Table 3. Material properties and dimensions of the four large arteries 



Element 
number 


Artery name 


Length 
(cm) 


Proximal area 
(cm^) 


Distal area 
(cm^) 


Elastic modulus 
(1 0^ dyne/cm) 


1 


Ascending aorta 


5.5 


6.605 


3.941 


4.0 


2 


Thoracic aorta 


18.5 


3.597 


2.835 


4.0 


3 


Abdominal aorta 


4.3 


2.378 


2.378 


4.0 


4 


Abdominal aorta 


9.6 


1.021 


1.021 


4.0 



total length of the four arteries was 37.9 cm, and the length of 
the thoracic aorta was about half (48.8%) of the total length, 
with a mean diameter of 2 cm. In the normal case, wall stiffness 
values for the arterial elements were identical to those de- 
scribed in the Model of Ozawa et al. (11). To couple the arterial 
network model with BLH, we regarded the left ventricular (LV) 
pressure as the boundary condition of the arterial network 
model and the aortic flow rate as the input condition of the 
BLH model. 

To investigate the effects of arterial wall stiffening on heart 
cellular mechanics and cardiovascular system dynamics, we 
performed a parametric study with respect to arterial wall stiff- 
ness, increasing arterial wall elasticity of the model by 1, 2, 3, 
and 4 times the normal value. However, this variation in arterial 
elasticity was assumed only in large arteries because the wall 
stif&iess of small arteries remains nearly unchanged (21). Here, 
PWV was calculated from the distance and traveling time of the 
pressure wave between the ascending aorta (1st element of the 
arterial network in Fig. 2) and the abdominal aorta (4th ele- 
ment of the arterial network in Fig. 2). 




Time (sec) 



Fig. 3. Numerical simulation results of sequential events of the integrative cell-sys- 
tem-arterial network model. Dotted lines, atrial events; solid lines, ventricular events. 
(A) Action potential. (B) Ca^+ concentration. (C) Cellular forces (solid line: active force 
generated by cross-bridges, thick dotted line: passive elastic force at the ventricular 
myocyte). (D) Cardiac volume. (E) Blood pressure in the left atrium and left ventricle, 
respectively (thick dotted line: aorta pressure). (F) Pressure pulse waves at three arte- 
rial positions (ascending aorta, common iliac artery, and anterior tibial artery). 

RESULTS 

To assess the feasibility of our cell-system-arterial network 
model, we plotted the sequential events from the cell level to 
that of the arterial pulse (Fig. 3) and reconstructed the shapes of 
the original action potential and Ca^^ transient. The action po- 
tential showed the characteristic spike notch dome shape of 
ventricular cells (Fig. 3A). Fig. 3B shows the transient increase 
in free calcium concentration within the cytoplasm. The force 
of muscle contraction develops rapidly in response to the Ca^^ 
transient in cells (Fig. 3C). Two components comprise this 
force: that generated by cross-bridges and elastic force. The LV 
active force showed a peak value of approximately 28 mN/mm^ 
(kPa). Fig. 3D shows the volume variation in the left atrium (LA) 
and LV relative to time. The ejection fraction of the left ventricle 
was nearly 0.58. The LA contracts prior to ventricular contrac- 
tion. The timeline of the pressure changes observed in the aor- 
ta, atrium, and ventricle are shown in Fig. 3E. Here, the heart 
rate was assumed to be 72 bpm (beats/min). Pressure pulse 
waves at three arterial positions (ascending aorta, common ili- 
ac artery, and anterior tibial artery) are plotted in Fig. 3F and 
show increased pulse pressures and a smoothed dicrotic notch 
in the distal portion of the arterial branches. 

To verify our model, the major hemodynamic variables gen- 
erated by the model are compared with standard values in Table 
4. The computed blood pressure is within the range considered 
to be physiologically normal within the general population. The 
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Fig. 4. Simulated results of propagation delay of arterial pressure waves between the ascending aorta and the abdominal aorta with respect to differences in arterial stiffness. 
(A) En, (B) 2En, (C) 4En. Here, En represents the normal stiffness of the large arteries. 



Table 4. Comparison between computational and normal parameter ranges 



Parameters 


Simulation 


Normal range 


References 


Pressure (mmHg), 








left ventricle 


117 


90-140 


(19) 


Systolic 


9 


4-12 


(19) 


End-diastolic 








Aorta 


116 


90-140 


(19) 


Systolic 


73 


60-90 


(19) 


Diastolic 








Elastic modulus of large 
arteries (1 0^ dyne/cm^) 


4 (normal case) 


2.1-6.1 


(21) 


PWV (m/s) 


4.8-8.3 


4.4-8.5 


(21) 



PWV, pulse wave velocity. 

stroke volume was approximately 73 mL. The cardiac output, 
which was obtained by multiplying the stroke volume with a 
heart rate of 72 bpm, was 5,256 mL/min. The computed value of 
the PWV in the aorta was within the range 4.4 to 8.5 m/sec, 
which is consistent with an experimental study in humans (22). 

Fig. 4 shows the variation in pressure waves within the aorta 
according to arterial stiffaess. Compared with the baseline val- 
ue of arterial elasticity (En = 4 x 10^ dyne/cm^), with increases 
in arterial stiffness, the time delay of the pressure wave between 
the ascending and abdominal aortas decreased, whereas the 
difference in peak aortic pressures increased. As shown in Fig. 
4A and B, increase of arterial stiffness by 2 times of normal val- 
ue increased peak aortic pressure by nearly 8% (from 112.5 
mmHg to 121.5 mmHg). However, the increased value is rela- 
tively small compared with in vivo experimental observation 
(23). Interestingly, in revers to the present cases, acute rise in 
blood pressure can increase in vivo arterial stiffaess demon- 
strated by Stewart et al. (24). The relation between arterial stiff- 
ness and cardiac output is shown in Table 5. 

At the cellular level, the cross-bridge elongation (CBEL) de- 
creased sharply during contraction below the steady-state value 
( = 5 nm) and then immediately reverted to the back to the 
steady-state value (Fig. 5A). Shortening of the CBEL below the 
steady-state value during contraction played a critical role in 
active force generation, as the calcium-bound troponin con- 



Table 5. Variations of cardiac output and pulse wave velocity according to the in- 
crease of arterial stiffness 



Elastic modulus 


Stroke 


Heart rate 


Cardiac output 


Pulse wave velocity 


(1 0^ dyne/cm) 


volume (ml) 


(beats/min) 


(mL/min) 


(m/s) 


4.0 


77.6 


72 


5,587 


4.8 


8.0 


75.9 


72 


5,465 


6.2 


16.0 


73.0 


72 


5,256 


8.3 



centration was high during this period (Fig. 5B). This shorten- 
ing of the elongation during contraction was accompanied by a 
decrease in the half sarcomere length (HSL) (Fig. 5C). The ac- 
tive force was generated by the product of the CBEL and the 
concentration of calcium-bound troponin (Fig. 5D). 

There was no change in the amount of troponin bound to the 
cross-bridge with increasing arterial thickness (solid, thick, and 
thin dotted lines in Fig. 5B), whereas CBEL increased (Fig. 5A). 
Accordingly, the HSL and LV peak pressure (Fig. 5E) increased 
with arterial stiffness. According to increased stiffness, it is re- 
markable that the changes of CBEL and half sarcomere length 
during ventricular relaxation were relatively small compared 
with those during ventricular contraction. 

The peak cellular stress observed with increases in arterial 
stiffness is shown in Fig. 6A. As shown in the Fig. 6, there was an 
initial steep increase in the peak cellular stress followed by a 
saturated pattern at higher arterial stiffness values. On the other 
hand, the variation in PWV was nearly linear with increasing 
arterial stiffness (Fig. 6B). Systolic blood pressures (SBP) in the 
ascending aorta and brachial artery also changed consistent 
with peak cellular stress (Fig. 6C). Interestingly, variation in aor- 
tic SBP was more prominent than that in the brachial arterial 
with respect to increasing arterial stiffness. Fig. 6D presents the 
relationship between PWV and peak cellular stress, showing a 
saturated pattern for higher values of PWV 

DISCUSSION 

A major objective of the present study was to integrate cellular 
processes with systemic and arterial hemodynamics by which 
propagation of the arterial pulse wave is related to heart me- 
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Fig. 5. Simulated results of the inotropic state of the 
left ventricular BLH in the cell-system model. (A) 
Cross-bridge elongation length (CBEL). (B) Concen- 
tration of cross-bridge-attached troponin ([TCa*]-^ 
[T*]). (C) Half sarcomere length. (D) Cross-bridge- 
generated contraction force (CBCF). (E) LV pressure 
according to arterial stiffness. Here, En represents 
the normal stiffness value of the large arteries. 
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chanics in the human circulatory system. Here, the arterial net- 
work model proposed by Ozawa et al. (11) could couple blood 
hemodynamics and pulse waves in detailed arteries with LV 
mechanics and systemic circulation. Our discussion of the 
computational results focuses on three main aspects: the plau- 
sibility of the integrative method of a cell-system- arterial net- 
work to simulate heart mechanics and cardiovascular hemody- 
namics, a parametric study focusing on arterial stiffness to de- 
lineate the relationship between arterial hemodynamics and 
cardiac cellular mechanics. 

In the cell-system- arterial network model, electrophysiologi- 
cal variables such as the action potential, calcium concentra- 
tion, and generated cellular forces in the cellular models were 
determined for atrial and ventricular cells (Fig. 3A-C). Compu- 
tational results for the cardiovascular system dynamics were 
consistent with published results of pressure and volume varia- 
tions in the LA and LV (Fig. 3D, E). Arterial pulse wave propaga- 
tion along the arterial network also showed a pattern similar to 
previous clinical observations (Fig. 3F, 4). In conclusion, all se- 
quential events from the cell to arterial pulse waves were well 
reproduced and consistent with the existing data, demonstrat- 
ing the plausibility of the present method as an integrative 
analysis tool for evaluating heart mechanics and arterial hemo- 
dynamics. 

Atherosclerosis is a symptom of aging that causes serious 
health problems. Specifically, it increases cardiac cellular stress, 
which eventually stimulates heart hypertrophy (23, 24). Roman 
et al. (25) investigated the effects of central blood pressure (BP) 
on left ventricular hypertrophy (LVH), and showed that systolic 
BP in the central artery (or the ascending aorta) is more impor- 
tant in stimulating LVH and inducing conformational changes 
and remodeling of the heart. With respect to medical terminol- 
ogy, central BP represents the BP in the large arteries, such as 
the ascending, descending, and abdominal aortas. Wang et al. 
(26) also reported that central systolic pressure is more critical 
than other blood pressures in predicting cardiovascular diseas- 
es. However, theoretical approaches to explain the detailed 
physiological mechanism underlying these effects have not 
been reported. Therefore, we proposed an integrative cell-sys- 
tem-arterial network model to delineate the relationships be- 
tween heart cellular stress, LVH, and arterial stiffaess. 

Using this model, we performed a parametric study to evalu- 
ate the effects of stiffened arteries on heart cellular mechanics 
by increasing the vascular wall stiffness of the largest arteries in 
the model. Cross-bridge-generated contraction force (CBCF) 
(Fig. 5D), LV peak pressure (Fig. 5E), and PWV (Fig. 6B) in- 
creased, whereas HSL decreased with increasing arterial stiff- 
ness. Peak cellular stress and SBP in the brachial and central ar- 
teries also increased, but were saturated at higher values of ar- 
terial stiffness (Fig. 6A). As shown in Fig. 6C, central arterial 
stiffness could more affect SBP, which is because the half of 



ejected blood is stored in the vessel and three quarters of the 
blood is stored in aorta. Therefore, if the central artery becomes 
stiff, the SBP could be more affected. Interestingly, changes in 
the peak cellular stress and SBP showed similar patterns (Fig. 
6A and C). In particular, SBP variation in the central artery (or 
ascending aorta) with respect to arterial stiffness was more 
prominent than that in the brachial artery) (Fig. 6C). As increas- 
es in exerted cellular stress indicate increased stimulation of 
LVH, this result clearly corresponds to the clinical observation 
that SBP in the central artery is more critical than other blood 
pressure measurements in predicting LVH. The peak cellular 
stress corresponding to a specific PWV can be inferred quanti- 
tatively from the graph shown in Fig. 6D. Clinically, cellular 
stress is not easy to observe in patients (27), while PWV can be 
measured easily. Therefore, this result provides another exam- 
ple of the utility of the present cell-system-arterial network 
model. 

This study developed a cell-system- arterial network coupled 
model of the cardiovascular system. Although results were rea- 
sonable, the study has several limitations. First, we did not con- 
sider other effects, such as autonomic nerve control or changes 
in the resistance of the peripheral vessels involved in the arteri- 
al stiffening process. Therefore, we elucidated only the mecha- 
nisms involved in arterial stiffness perturbation, cellular con- 
tractility, and vascular hemodynamics. Second, the heart has a 
complex muscle layer with fiber anisotropy and heterogeneous 
muscle thickness. In addition, different parts of the ventricles 
are activated at different times. However, in this study, the ven- 
tricle was assumed to be a Laplace heart to reduce the com- 
plexity of the coupled cell- circulation hemodynamics method. 
Besides, in this study we didn't consider long-term compensa- 
tory mechanism of blood pressure through renal excretory sys- 
tem. Nevertheless, these limitations were not expected to mark- 
edly alter the main findings of this study. 

Despite its limitations, our model can simulate the sequen- 
tial events of the cardiovascular system from cells to the sys- 
temic circulation. Consequently, the model can be applied to 
perioperative intensive care in heart, lung and liver surgery 
which requires accurate evaluation of the pulmonary preload 
and systemic afterload. Along with clinical data, a computer 
model of the cardiovascular system can be a useful tool allow- 
ing clinicians to estimate critical information about patients' 
cardiovascular system, such as blood pressure, and the arterial 
pulse wave. In addition, clinicians can use the model to predict 
the cardiovascular changes induced in patients by treatment in 
intensive care units. For example, Kashif et al. (28) showed that 
a simulation method combined with clinical measurements 
can provide information about patients in intensive care. We 
expect that our model could play a similar role in estimating 
patient-specific data for grave patients. In the future, this will be 
delineated by combining the model with clinical data. 
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